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“ iAn  experimental  program  was  conducted  to  investigate  the 
interaction  effects  which  occur  between  the  nozzle  exhaust  flow  and 
the  external  fiow  field  associated  with  isolated  nozzle  afterbody 
configurations  at  transonic  Mach  numbers.  Pressure  data  were 
obtained  for  three  afterbody  geometries  with  boattail  angles  of 
10,  15,  and  25  deg  at  Mach  numbers  from  0.6  to  1.5  at  zero  angles 
of  attack  and  sideslip.  Cold  (high-pressure  air)  and  hot 
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(air/ethylene  combustion)  jet  test  techniques  were  used  to  simulate 
and  duplicate,  respectively,  the  nozzle  exhaust  flow  for  a sonic 
jet  installation.  Nozzle  exhaust  temperature  was  varied  from  540 
to  approximately  2,900^R.  The  most  significant  results  pertain  to 
those  effects  on  boattail  pressure  drag  caused  by  exhaust  plume 
temperature  and  flow  asymmetry  (model  support  strut  induced) . The 
differences  obtained  in  boattail  pressure  drag  between  the  cold 
jet  simulation  and  hot  jet  duplication  results  were  significant  at 
nozzle  pressure  ratios  representative  for  turbofan  and  turbojet 
engines  at  subsonic  Mach  numbers.  Adjusting  the  cold  jet  nozzle 
pressure  ratio  to  correct  for  changes  in  the  exhaust  plume  specific 
heat  ratio  with  temperature  did  not  account  for  the  differences 
observed.  Flow  asymmetry  effects  were  Mach  number  and  nozzle 
pressure  ratio  dependent  and  increased  in  severity  as  the  boattail 
angle  was  increased. 
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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEdC),  Air  Force  Systems  Command  (AFSC),  under  Program  Element  65807F. 
The  results  of  the  tests  were  obtained  by  ARO,  Inc.  (a  subsidiary  of  Sverdrup  & Parcel 
and  Associates,  Inc.),  contract  operator  of  AEDC,  AFSC,  Arnold  Air  Force  Station, 
Tennessee,  under  ARO  Project  Number  PA411^  in  support  of  Research  Project  Number 
PF430  and  under  ARO  Project  Number  P41T-80A  in  support  of  Research  Project  Numbers 
P32P^1 1 B and  P32P-AOA.  The  authors  of  this  report  were  L.  L.  Galigher,  S.  F.  Yaros, 
and  R.  C.  Bauer,  ARO,  Inc.  The  manuscript  (ARO  Control  No.  ARO-PWT-TR-76-28)  was 
submitted  for  publication  on  March  9,  1976. 
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1.0  INTRODUCTION 

The  technical  challenge  associated  with  exhaust  system  integration  is  formidable  in 
view  of  the  complex  interactions  which  occur  between  the  exhaust  plume,  the  local  free 
stream,  and  the  afterbody  surfaces.  The  flow  interactions  between  the  exhaust  plume  and 
the  free  stream  can  be  categorized  as  either  inviscid  or  viscous  in  nature.  Inviscid 
interactions  result  if  the  jet  exhaust  is  not  fully  expanded  and  the  local  external  flow 
does  not  match  the  jet  exhaust  Mach  number  and  flow  angle  - a plume  displacement 
effect.  Viscous  interactions  result  from  mixing  at  the  interface  of  the  local  external  flow 
and  the  jet  exhaust.  For  simulation  of  the  flow  field  about  the  nozzle  afterbody,  wind 
tunnel  test  programs  normally  use  high-pressure  air  to  simulate  nozzle  exhaust  flow  - not 
only  for  economic  reasons  but  also  to  simplify  testing  requirements.  The  cold  jet  technique 
can  be  used  to  match  either  the  initial  inclination  angle  or  the  maximum  diameter  of 
the  inviscid  plume  boundary  of  the  hot  jet.  However,  viscous  interaction  effects  for  a 
cold  jet  differ  from  those  for  a hot  jet  since  mixing  between  the  nozzle  exhaust  and 
local  external  flow  depends  largely  on  temperature  and  velocity  gradients. 

The  purpose  of  the  experimental  program  reported  herein  was  to  investigate  the 
interaction  effects  which  occur  between  the  nozzle  exhaust  flow  and  the  external  flow 
field  associated  with  isolated  nozzle/afterbody  (NAB)  installations  in  the  transonic  Mach 
number  range.  Configuration  variables  included  nozzle  boattail  geometry  and  internal  nozzle 
geometry.  A cold  jet  and  hot  jet  test  technique  was  used  to  simulate  and  duplicate, 
respectively,  the  nozzle  exhaust  flow.  High-pressure  air  served  as  the  simulation  fluid,  and 
an  air/ethylene  combustor  provided  the  means  for  duplicating  the  exhaust  flow  for  the 
hot  jet  installation.  The  experimental  program  concerning  the  NAB  geometry  effects  was 
initiated  in  response  to  a request  by  the  Advisory  Group  for  Aeronautical  Research  and 
Development  (AGARD)  of  NATO  for  member  nations  to  participate  in  a research  project 
on  improved  nozzle  testing  techniques  in  transonic  flow.  The  overall  objective  of  the 
AGARD  program  is  to  parametrically  develop  the  influences  of  exhaust  nozzle  flows  on 
the  external  afterbody  flow  and,  conversely,  show  the  effect  of  external  afterbody  flow 
on  nozzle  internal  flow  characteristics.  The  hot  jet  duplication  experiments  were  conducted 
by  AEDC  as  part  of  an  effort  to  investigate  hot  versus  cold  jet  effects  on  boattails  of 
different  geometries. 


2.0  APPARATUS 


2.1  TEST  FACILITY 

The  transonic  Propulsion  Wind  Tunnel  (16T)  is  a continuous  flow,  closed-circuit  wind 
tunnel  capable  of  operating  within  a Mach  number  range  from  0.2  to  1 .6.  The  tunnel 
can  be  operated  within  a stagnation  pressure  range  from  120  to  4,000  psfa  depending 
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on  the  Mach  number.  The  stagnation  temperature  can  be  varied  from  an  average  minimum 
of  about  80  to  a maximum  of  160°F  depending  on  the  cooling  water  temperature.  The 
specific  humidity  of  the  air  in  the  tunnel  is  controlled  by  exchanging  tunnel  air  for 
conditioned  makeup  air  from  an  atmospheric  drier. 

2.2  EXPERIMENTAL  HARDWARE  AND  MODEL  SUPPORT  SYSTEM 

An  existing  axisymmetric  cone-cylinder  model  was  modified  to  satisfy  requirements 
for  both  cold  and  hot  exhaust  testing  with  three  external  afterbody  configurations  having 
boattail  angles  of  10,  15,  and  25  deg.  The  model  had  a length  of  approximately  147 
in.,  a maximum  diameter  of  9.86  in.,  and  a 14-deg,  half-angle  conical  nose.  A 
boundary-layer  trip  consisting  of  Q.055-in.-diam  steel  spheres  spot  welded  to  a trip  ring 
at  a circumferential  spacing  of  four  sphere  diameters  was  located  12  in.  aft  of  the  cone 
vertex.  The  model  was  mounted  on  a tapered  strut  with  an  average  aft  sweep  angle  of 
35  deg.  The  strut  thickness-to-chord  ratio  varied  from  0.053  at  the  model  to  0.088  at 
the  tunnel  floor.  The  maximum  cross-sectional  area  of  the  model/strut  arrangement  was 
equivalent  to  0.88  percent  of  the  wind  tunnel  test  section  cross-sectional  area.  A 
dimensioned  sketch  of  the  model  is  presented  in  Fig.  1,  and  the  cross-sectional  area 
distribution  of  the  model/strut  arrangement  is  shown  in  Fig.  2.  A dimensioned  sketch 
of  each  of  the  three  nozzle  boattail  configurations  (10-,  15-,  and  25-deg  boattail  angle) 
is  presented  in  Fig.  3.  An  installation  photograph  of  the  model  installed  in  Tunnel  16T 
is  presented  in  Fig.  4. 

2.2.1  Cold  Flow  Nozzle  Assembly 

High-pressure  air  was  used  to  simulate  the  jet  exhaust  for  the  cold  flow  portion  of 
the  experiments.  The  high-pressure  air  was  ducted  through  the  model  support  strut  to 
a plenum  in  the  model  forebody.  The  air  continued  aft  through  a flow-conditioning  section 
prior  to  entering  the  convergent  nozzle  assembly.  The  nozzle  assembly,  which  had  a 
contraction  ratio  of  3.28  and  an  exit  diameter  of  3.982  in.,  is  sketched  in  Fig.  5. 

2.2.2  Hot  Flow  Nozzle  Assembly 

An  ethylene/air  combustor  was  used  to  provide  hot  jet  exhaust  flow  duplication. 
Ethylene  (C2  H4  ) is  a gaseous  hydrocarbon  fuel  which,  when  burned  in  air,  produces  exhaust 
products  that  very  closely  duplicate  the  exhaust  products  of  JP-4  burned  in  air.  The 
water-cooled  combustor  assembly  was  fabricated  from  copper  and  had  an  inside  diameter 
of  6.15  in.  and  a wall  thickness  of  0.125  in.  The  combustor  terminated  with  a convergent 
nozzle  with  an  exit  diameter  of  3.228  in.  and  a contraction  ratio  of  3.63.  The  flame 
holder  consisted  of  four  doughnut  rings  interconnected  with  fuel  flow  passages,  but  only 
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the  forward  doughnut  ring  served  as  the  ethylene  injection  ring.  Ignition  was  accomplished 
by  injecting  a small  quantity  of  tri-ethyl  borane  (TEB,  a pyrophoric  fuel)  into  a retainer 
mounted  or.  the  flame  holder.  Flow  mixers  were  placed  at  the  rear  of  the  flame  holder 
to  mix  the  hot  flow  from  the  combustion  zone  with  the  air  around  the  outside  of  the 
flame  holder.  A sketch  of  the  hot  flow  nozzle  assembly  is  shown  in  Fig.  6.  An  installation 
photograph  of  the  hot  flow  model  is  presented  in  Fig.  7.  The  nozzle  base  area  depicted 
in  the  preceding  photograph  was  seven  times  larger  than  that  for  the  cold  flow  model 
(4.968  in.2  compared  to  0.709  in.2).  This  difference  in  nozzle  base  area  resulted  from 
imposing  water  cooling  requirements  on  the  hot  flow  nozzle  assembly. 

2.3  INSTRUMENTATION 

External  pressure  instrumentation  was  common  for  both  the  cold  and  hot  flow  model 
configurations;  only  the  instrumentation  within  the  nozzle  flow  chambers  was  different. 
Rows  of  static  pressure  orifices  were  located  along  the  model  centerbody  and  on  the 
nozzle  boattails  at  various  circumferential  stations  as  shown  in  Table  1.  Pressure  and 
temperature  instrumentation  in  the  nozzle  flow  chambers  of  the  cold  and  hot  flow  model 
configurations  are  presented  in  Tables  2 and  3,  respectively.  Air  and  ethylene  flow  rates 
were  determined  from  pressure,  temperature,  and  area  measurements  in  critical  flow  venturi 
metering  sections  external  to  the  models. 

3.0  PROCEDURE 

3.1  EXPERIMENTAL  CONDITIONS  AND  TECHNIQUE 

The  data  presented  in  this  report  were  obtained  at  nominal  free-stream  Mach  numbers 
of  0.6,  0.8,  0.9,  0.95,  1.1,  and  1.5  at  a free-stream  Reynolds  number  of  2.5  x 106  per 
foot.  Data  on  the  influence  of  Reynolds  number  for  the  range  from  1 to  50  x 106  per 
foot  are  reported  in  Ref.  1 . Angle  of  attack  was  zero  at  all  test  conditions.  Nozzle  total 
pressure-to-free-stream  static  pressure  ratio  (NPR)  was  varied  from  jet-off  conditions  to 
a maximum  of  22  depending  on  the  free-stream  Mach  number.  Nozzle  pressure  ratio  surveys 
for  the  hot  flow  model  were  conducted  primarily  at  constant  fuel/air  ratios  of  0,  0.015, 
0.025,  and  0.040.  The  10-,  15-,  and  25-deg  boattail  configurations  were  tested  on  the 
cold  flow  model,  but  only  the  1 5-  and  25-deg  boattail  configurations  were  tested  on  the 
hot  flow  model.  Schlieren  photographs  were  obtained  at  selected  nozzle  pressure  ratios 
for  both  the  cold  and  hot  flow  models. 

The  data  acquisition  procedure  consisted  of  (1)  setting  desired  free-stream  Mach 
number,  (2)  obtaining  jet-off  data,  (3)  setting  desired  fuel/air  ratio,  if  applicable,  and 
varying  nozzle  total  pressure  ratio,  and  (4)  obtaining  schlieren  photographs  as  required. 
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32  DATA  REDUCTION 


The  primary  performance  parameter  presented  in  this  report  is  the  nozzle  boattail 
pressure  drag  coefficient  obtained  from  integration  of  the  experimental  pressure 
distributions.  The  data  reduction  technique  employed  was  based  on  using  the  measurements 
of  each  individual  row  of  pressure  orifices  located  at  circumferential  stations  of  0,  180, 
and  270  deg  and  assuming  that  each  of  these  pressure  distributions  was  imposed  on  all 
of  the  boattail  projected  area.  Thus,  pressure  drag  coefficients  based  on  the  top,  side, 
and  bottom  row  of  pressure  orifices  were  calculated  as  follows: 


CDPX 


(Note:  X = top,  iide,  or  bottom) 


where  "A"  is  the  reference  area  based  on  the  maximum  diameter  of  the  model,  "CPi" 
is  the  pressure  coefficient  for  orifices  identified  in  Table  1,  and  "Ri"  is  the  radius  of 
the  boattail  surface  relative  to  the  nozzle  centerline.  The  drag  coefficient  values,  as 
calculated,  pertain  to  the  pressure  drag  on  the  nozzle  boattail  projected  area  aft  of  model 
station  130.471  in.  (see  Fig.  1). 


For  the  cold  flow  nozzle  assembly,  the  nozzle  total  pressure,  PTJ,  was  calculated 
from  one-dimensional  relationships  using  measured  nozzle  mass  flow  and  nozzle  flow  duct 
static  pressure,  total  temperature,  and  cross-sectional  area.  Nozzle  discharge  coefficient, 
CDND,  was  calculated  by  the  following  equation: 


CDND  = 1.8803 


(PTJ)(AT)(TTJ)'0,5 


For  the  hot  flow  nozzle  assembly,  static  pressure  orifices  located  within  the  flow 
duct  and  just  forward  of  the  start  of  the  throat  contraction  were  used  to  calculate  the 
nozzle  total  pressure,  PTJ,  in  the  following  manner: 


PTJ  = P, 
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where  (PT/P)i  is  the  subsonic  isentropic  pressure  ratio  based  on  the  nozzle  contraction 
ratio  of  3.63  and  the  ratio  of  specific  heats,  y , of  the  exhaust  products  based  on  100-percent 
combustion  efficiency.  Since  PT/P  is  an  extremely  weak  function  of  y (PT/P  decreases 
approximately  0.15  percent  for  a change  in  y from  1.40  to  1.26),  use  of  y based  on 
100-percent  combustion  efficiency  introduces  no  significant  error  in  the  calculation  of 
PTJ.  Exhaust  total  temperature  was  calculated  from  mass  flow  relationship  in  the  following 
manner: 


i 

’Kv(PTJ)(AT)(Cr)NI))T 

TTJ  . : : 

L M,  + Mf  J 


where  Ky  is  a function  of  y,  AT  is  the  measured  throat  area,  CDND  is  the  nozzle  discharge 
coefficient  (0.9902)  determined  from  cold  flow  tests  of  the  hot  flow  nozzle,  and  Ma  and  Mf 
are  air  and  ethylene  mass  flow  rates,  respectively. 


3.3  PRECISION  OF  MEASUREMENTS 


The  uncertainty  of  the  major  test  parameters  is  tabulated  below.  The  uncertainty 
in  these  parameters  includes  the  inaccuracies  in  the  tunnel  reference  systems,  the  recording 
systems,  and  the  pressure  transducers  themselves.  Uncertainty  is  defined  as 


r « *(B  «-  2S) 


where  "B"  is  the  bias  and  "S"  is  an  estimate  of  the  standard  deviation,  both  values  being 
determined  from  the  Taylor  series  propagation-of-error  method. 

UNCERTAINTY 

Mach  Number  . 


Parameter 

0.6 

0.8 

0.9 

0.95 

1.1 

1.5 

M„ 

0.00320 

0.00300 

0.00300 

0.00290 

0.00300 

0.00360 

CDPX 

0.00056 

0.00046 

0.00043 

0.00038 

0.00035 

0.00031 

CP 

0.00460 

0.00330 

0.00290 

0.00270 

0.00240 

0.00200 

NPR 

0.00380 

0.00510 

0.00600 

0.00630 

0.00760 

0.01290 

4.0  EXPERIMENTAL  RESULTS 


The  experimental  results  are  presented  in  ‘he  following  order: 

4.1  Flow  Duct  and  Nozzle  Entrance  Characteristics 

4.2  Nozzle  Pressure  Ratio  Effects 

4.3  Mach  Number  Effects 
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4.4  Flow  Asymmetry  Effects 

4.5  Exhaust  Plume  Temperature  Effects 

4.6  Nozzle  Base  Area  Effects 

The  cold  flow  model  data  are  presented  to  illustrate  the  effects  of  nozzle  pressure  ratio, 
Mach  number,  and  flow  asymmetry  on  boattail  pressure  distribution  and  drag  coefficient 
since  experiments  on  all  three  boattail  configurations  were  conducted  during  the  cold  flow 
phase.  As  stated  previously,  experiments  were  conducted  only  on  the  15-  and  25-deg 
boattail  configurations  with  the  hot  flow  model.  The  exhaust  temperature  effects  were 
determined  by  comparing  data  from  the  hot  flow  model  operated  with  and  without  the 
combustor  in  operation. 

The  boattail  pressure  drag  coefficient,  as  presented  in  this  report,  was 
nondimensionalized  with  the  model  maximum  cross-sectional  area  and  should  not  be 
construed  as  being  equal  in  magnitude  to  aircraft  drag  counts  (one  drag  count  = 0.0001 
in  drag  coefficient),  which  are  based  on  wing  area.  (Generally  an  aircraft  drag  count  for 
an  aircraft  similar  to  a lightweight  fighter  is  equivalent  to  approximately  15  "body"  drag 
counts.) 

4.1  FLOW  DUCT  AND  NOZZLE  ENTRANCE  CHARACTERISTICS 

Nozzle  total  pressure  for  the  cold  flow  model  was  determined  by  two  methods:  one 
from  one-dimensional  relationships  using  the  measured  mass  flow  and  flow  duct  static 
pressure,  total  temperature,  and  cross-sectional  area,  and  the  other  from  a ten-probe, 
area-weighted,  total  pressure  rake  installed  in  the  flow  duct.  The  flow  duct  total  pressure 
profiles  shown  in  Fig.  8 were  representative  of  those  profiles  established  for  various  flow 
duct  pressure  levels.  The  total  pressure  distortion  was  approximately  5.4  percent  and  was 
independent  of  flow  duct  pressure  level.  The  difference  between  the  calculated  total 
pressure  value  and  the  measured  total  pressure  value  obtained  from  the  rake  was  less  than 
1 percent.  The  data  presented  in  Fig.  9 show  the  excellent  agreement  obtained  between 
the  two  methods  used  to  obtain  flow  duct  total  pressure.  For  this  test,  the  calculated 
flow  duct  total  pressure  was  considered  the  primary  measurement,  and  the  average  total 
pressure,  as  determined  from  the  pressure  rake,  was  considered  as  the  backup  measurement. 

The  discharge  coefficient  (CDND)  for  the  cold  flow  model  is  presented  in  Fig.  10 
as  a function  of  throat  Reynolds  number  (RET).  The  discharge  coefficient  generally  varied 
from  approximately  0.994  at  a throat  Reynolds  number  of  I x 106  to  approximately 
0.998  at  a throat  Reynolds  number  of  5 x 106.  For  the  hot  flow  model,  the  nozzle 

discharge  coefficient  had  an  average  value  of  0.99. 

i 

j 

1 The  calculated  values  of  exhaust  total  temperature  for  the  hot  flow  model  are 

presented  in  Fig.  1 1 as  a function  of  flow  duct  total  pressure.  These  data  indicate  that 
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combustion  efficiency  is  pressure  level  dependent,  especially  at  a fuel/air  ratio  of  0.04. 

The  average  exhaust  total  temperatures  at  fuel/air  ratios  of  0.015,  0.025,  and  0.04  were 
approximately  1,600,  2,100,  and  2,900°R,  respectively.  ' 

4.2  NOZZLE  PRESSURE  RATIO  EFFECTS 

t 

» 

The  basic  data  showing  the  effects  of  nozzle  pressure  ratio  on  the  pressure  drag  and  j 
pressure  distributions  for  the  10-,  15-,  and  25-deg  boattail  configurations  are  presented 
in  Figs.  12  through  15.  The  data  presented  in  these  figures  represent  the  data  obtained 
from  the  top  row  of  pressure  orifices.  In  general,  the  drag  trend  with  nozzle  pressure 
ratio  (NPR),  as  well  as  the  drag  level,  for  the  25-deg  boattail  configuration  differed  i 

significantly  from  that  for  the  10-  and  15-deg  boattail  configurations  at  subsonic  Mach  i 

numbers;  however,  as  the  Mach  number  increased  supersonically,  all  three  boattail 
configurations  exhibited  similar  drag  trends,  and  the  differences  in  drag  level  decreased  i 
noticcahly.  < 

The  pressure  distribution  data  presented  for  the  10-  and  15-deg  boattail  configurations  | 
are  characteristic  of  these  distributions  usually  encountered  with  small  angle  boattails  at 
subsonic  Mach  numbers  - flow  expansion  around  the  boattail  shoulder  followed  by  a 
well-behaved  flow  recompression  to  a pressure  level  above  free-stream  ambient.  The  pressure 
distributions  for  the  25-deg  boattail  configuration  at  similar  Mach  numbers  differed 
radically  from  those  obtained  for  the  10-  and  15-deg  configurations  in  that  significantly  j 
more  flow  expansion  occurred  and  was  followed  by  an  abrupt  flow  recompression  to  a 
pressure  level  approacliing  the  ambient.  It  is  difficult  to  discern  from  the  schlieren 
photographs  presented  in  Fig.  1 6 if  the  very  adverse  pressure  gradient  (AP/AX  > 0)  imposed 
on  the  25-deg  boattail  configuration  at  subsonic  Mach  numbers  resulted  in  flow  separation;  , 

however,  it  is  suspected  that  the  adverse  pressure  gradient  imposed  on  the  25-deg  boattail  j 

configuration  was  of  sufficient  magnitude  to  cause  separation  of  the  turbulent  boundary  j 

layer.  j 

At  supersonic  Mach  numbers,  only  the  pressure  distribution  data  for  the  10-deg  ! 
boattail  configuration  show  flow  recompression  above  free-stream  static  pressure,  and 
compared  with  the  subsonic  data,  the  flow  recompression  was  much  less.  The  schlieren 
photographs  obtained  at  Mach  number  1.5  show  quite  clearly  the  shock  wave  structure 
in  the  flow  region  near  the  10-  and  15-deg  boattail  surfaces.  Flow  separation  is  visible  j 

in  the  schlieren  photograph  for  the  25-deg  boattail  configuration  at  Mach  number  1.5.  i 

4.3  MACH  NUMBER  EFFECTS  ' 

The  data  presented  to  show  Mach  number  effects  on  boattail  drag  coefficient  and 
pressure  distribution  are  for  scheduled  nozzle  pressure  ratio  values  of  3 at  all  subsonic 
Mach  numbers  and  5 and  7 for  M„  = 1.1  and  1.5,  respectively.  The  most  obvious  aspect 
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of  the  drag  coefficient  data  presented  in  Fig.  1 7 was  the  onset  of  drag  rise  characteristic 
exhibited  for  each  of  the  boattail  configurations.  As  boattail  angle  was  increased  from 
10  to  25  deg,  the  Mach  number  at  which  the  drag  rise  occurred  decreased.  The  drag 
rise  occurred  at  = 0.9  for  the  10-deg  boattail  configuration,  at  = 0.8  for  the  15-deg 
boattail  configuration,  and  at  a Mach  number  equal  to  or  less  than  0.6  for  the  25-deg 
boattail  configuration. 

Before  the  onset  of  drag  rise  for  the  10-  and  15-deg  boattail  configurations,  the  drag 
coefficient  levels  were  essentially  at  a constant  value.  The  pressure  distributions  for  the 
10-  and  15-deg  boattail  configurations,  as  shown  in  Figs.  18  and  19,  at  Mach  numbers 
less  than  that  corresponding  to  the  onset  of  drag  rise  show  that  although  flow  expansion 
around  the  shoulder  of  the  boattails  increased  as  Mach  number  was  increased,  flow 
recompression  also  increased.  The  net  result  was  that  no  change  in  pressure  drag  occurred. 

4.4  FLOW  ASYMMETRY  EFFECTS 

Flow  asymmetry  effects  were  determined  from  independent  calculations  of  boattail 
pressure  drag  based  on  the  pressure  distribution  as  measured  from  the  (1)  top  row  of 
pressure  orifices,  (2)  side  row  of  pressure  orifices,  and  3)  bottom  row  of  pressure  orifices. 
The  bottom  row  of  pressure  orifices  was  in  the  flow  field  of  the  model  support  strut 
wake.  Any  differences  obtained  in  the  drag  values  calculated  in  the  above  manner  were 
the  result  of  flow  asymmetry.  Nozzle  pressure  ratio  effects  on  flow  asymmetry  are  shown 
in  Figs.  20  through  22  for  the  10-,  15-,  and  25-deg  boattail  configurations,  respectively. 
The  drag  coefficient  data  presented  in  Fig.  23  summarize  flow  asymmetry  results  as  a 
function  of  Mach  number  (at  the  NPR  schedule  defined  in  Section  4.3).  Pressure 
distribution  data  for  each  of  the  boattail  configurations  are  presented  in  Figs.  24  through 
26. 


In  general,  flow  asymmetry  effects  were  Mach  number  and  NPR  dependent;  this  was 
especially  true  for  the  25-deg  boattail  configuration,  less  so  for  the  15-deg  boattail 
configuration,  and  practically  not  at  all  for  the  10-deg  boattail  configuration.  For  a given 
NPR,  flow  asymmetry  effects  increased  as  boattail  angle  was  increased.  The  drag  coefficient, 
as  determined  from  pressure  integration  of  the  bottom  row  of  pressure  orifices  for  the 
25-deg  boattail  configuration,  was  significantly  larger  than  that  obtained  for  the- top  row 
of  pressure  orifices.  For  example,  at  M^  = 0.95  and  for  NPR  = 3 the  drag  coefficient 
as  determined  from  the  bottom  row  of  pressure  orifices  was  approximately  16  percent 
larger  than  that  for  the  top  row  of  pressure  orifices.  This  16-percent  difference,  when 
translated  into  aircraft  drag  counts,  is  equivalent  to  approximately  7 percent  of  the  aircraft 
drag  based  on  a total  drag  of  250  counts.  Similar  analysis  for  the  15-deg  boattail 
configuration  at  M„  = 0.95  would  result  in  4 percent  of  the  aircraft  drag.  For  the  10-deg 
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boattail  configuration,  the  analysis  would  result  in  zero  percent  of  the  aircraft  drag  since 
the  pressure  drag  coefficients  based  on  the  top  and  bottom  row  of  pressure  orifices  were 
equal.  The  pressure  distribution  data  for  the  25-deg  boattail  configuration  show  a significant 
difference  between  the  top  and  bottom  rows  of  pressure  orifices  - on  the  order  of  0.04 
in  pressure  coefficient  aft  of  the  boattail  shoulder.  It  should  be  noted  that  a 0.01  change 
in  the  average  pressure  coefficient  for  the  boattails  is  equivalent  to  approximately  0.0085 
in  boattail  pressure  drag  coefficient. 

4.5  EXHAUST  PLUME  TEMPERATURE  EFFECTS 

The  hot  flow  model  data  were  obtained  at  fuel/air  ratios  of  0,  0.015,  0.025,  and 
0.04,  which  corresponded  to  exhaust  gas  temperatures  of  approximately  540,  1,600,  2,100, 
and  2,900°R,  respectively.  Over  this  temperatuie  range,  the  ratio  of  specific  heats  (y)  varied 
from  1.4  to  approximately  1.28.  The  exhaust  plume  temperature  effects  on  the  boattail 
drag  for  the  15-  and  25-deg  boattail  configurations  are  shown  in  Figs.  27  and  28, 
respectively.  These  data  illustrate  the  fact  that  "cold  jet"  test  results  must  first  be  adjusted 
for  exhaust  temperature  effects  before  they  car,  be  used  to  estimate  nozzle  boattail 
performance  for  a "live"  engine  installation.  A procedure  normally  used  to  correct  "cold 
jet"  test  results  is  based  on  the  premise  that  the  jet  plume  characteristics  immediately 
downstream  of  the  nozzle  exit  are  functions  primarily  of  nozzle  geometry,  nozzle  total 
pressure  ratio,  and  the  ratio  of  specific  heats.  The  effect  of  increasing  the  nozzle  total 
pressure  ratio  is  to  increase  the  initial  inclination  angle  of  the  jet  and  to  increase  the 
maximum  plume  diameter.  The  effect  of  increasing  y of  the  jet  is  to  decrease  the  initial 
inclination  angle  of  the  jet  and  to  decrease  the  maximum  plume  diameter.  Simulation 
of  a "hot  jet"  (y  < 1.4)  plume  operating  at  a given  nozzle  pressure  ratio  using  cold 
air  requires  that  the  nozzle  pressure  ratio  be  increased  to  offset  the  effect  of  y.  For  example, 
to  simulate  the  initial  inclination  angle  of  a "hot  jet"  (with  sonic  exit  conditions)  operating 
at  a nozzle  total  pressure  ratio  of  4.5  and  with  a value  of  y of  1.286,  the  "cold  jet" 
must  be  operated  at  a nozzle  total  pressure  ratio  of  5. 

If  the  procedure  as  outlined  above  and  described  in  Ref.  2 is  correct,  then  the  data 
presented  in  Figs.  27  and  28  for  the  1 5-  and  25-deg  boattail  configurations  should  collapse 
into  one  curve  when  plotted  as  a function  of  initial  jet  inclination  angle  (called  "y 
correction"  herein).  As  shown  in  Figs.  29  and  30,  the  "y  correction"  was  dependent  on 
Mach  number,  nozzle  pressure  ratio,  and  boattail  configuration.  In  general,  the  data 
collapsed  into  one  curve  only  at  = 1.5  over  the  range  of  nozzle  pressure  ratio  tested. 
At  Mach  numbers  of  less  than  1.5,  the  "y  correction"  tended  to  collapse  the  data  only 
at  the  higher  nozzle  pressure  ratio  values.  For  the  25-deg  boattail  configuration  at  subsonic 
Mach  numbers,  the  data  do  not  tend  to  collapse  at  any  nozzle  pressure  ratios  tested. 
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It  should  be  noted  that  application  of  the  "7  correction"  to  "cold  jet"  test  results 
at  nozzle  pressure  ratios  representative  of  turbofan  and  turbojet  engines  can  result  in 
significant  error  if  the  "corrected  cold  jet"  drag  values  are  used  to  estimate  boattail 
performance  for  "live"  engine  installations.  For  example,  at  = 0.8  and  NPR  - 3,  use 
of  the  "7  correction"  technique  to  correct  15-deg  boattail  "cold  jet"  test  results  for 
application  to  a "live"  engine  operating  with  a plume  temperature  of  2,900°R  would  result 
in  a 45-percent  error  (increase)  in  boattail  drag  coefficient.  When  translated  into  aircraft 
drag  counts,  this  increase  is  equivalent  to  approximately  3 percent  of  the  total  aircraft 
drag  (based  on  a total  aircraft  drag  of  250  drag  counts).  Additional  results  of  correlating 
plume  temperature  effects  are  reported  in  Ref.  3. 

4.6  NOZZLE  BASE  AREA  EFFECTS 


Water  cooling  requirements  for  the  hot  flow  nozzle  assembly  resulted  in  a nozzle 
base  area  equivalent  to  6.5  percent  of  the  model  maximum  cross-sectional  area,  as  compared 
to  0.93  percent  for  the  cold  flow  nozzle  assembly.  The  data  presented  in  Fig.  31  for 
the  15-  and  25-deg  boattail  configurations  show  that  nozzle  base  area  effects  on  boattail 
drag  were  significant  and  dependent  on  nozzle  pressure  ratio  and  Mach  number.  For  the 
15-deg  boattail  configuration,  the  drag  coefficient  as  measured  on  the  hot  flow  model 
was  either  equal  to  or  significantly  larger  than  that  measured  on  the  cold  flow  model 
through  the  Mach  number  range  from  0.6  to  1.5.  The  25-deg  boattail  configuration 
exhibited  similar  characteristics  except  at  Mw  = 0.95,  where  the  drag  as  measured  on  the 
cold  flow  model  at  nozzle  pressure  ratios  of  less  than  6 was  significantly  larger  than  that 
measured  on  the  hot  flow  model.  The  following  table  summarizes  base  area  effects  on 
the  drag  coefficients  for  the  15-  and  25-deg  boattail  configurations  at  the  nozzle  pressure 
ratio  schedule  defined  previously  in  Section  4.3. 


Percent 


(,'DI’T  (Hot  Flow  Model)  - COPT  (Cold  Flow  Model) 
COPT  (Cold  Flow  Model) 


x 100 


15-deg  Boattail.  25-deg  Boattail, 

M^  NPR  percent  percent 
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3 
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5.0  THEORETICAL  ANALYSIS 

5.1  PHYSICAL  STATEMENT  OF  PROBLEM 

One  of  the  main  concerns  of  this  analysis  was  the  accurate  description  of  the  physical  ; 

phenomena  involved  throughout  the  complex  flow  field.  The  analytical  model  chosen  was 
a compromise  between  the  requirements  of  the  physical  reality  and  the  capabilities  of  ; 

the  available  analytical  methods.  Since  much  was  discovered  concerning  both  facets  of  » 

this  problen  during  the  course  of  the  analyses,  the  analytical  model  was  constantly  revised 
and  is,  indeed,  still  being  modified.  In  general  the  model  accurately  represented  a large 
portion  of  the  flt«w  and  allowed  fairly  accurate  predictions  of  the  pressure  distribution  s 

over  the  nozzle  boattail.  The  following  discussions  of  the  actual  flow  characteristics  and 
the  analytical  model  should  indicate  the  validity  of  the  model  and,  to  a certain  extent,  ; 

the  areas  of  disagreement  between  theory  and  experiment. 

5.1.1  Characteristics  of  the  Actual  Flow 

The  configurations  that  were  analyzed  were  assumed  to  be  nonseparated  with  thin  j 
boundary  layers  (relative  to  the  reference  maximum  radius  of  the  cone-cylinder  model).  j 

Although  this  was  the  case  for  most  of  the  experiments  conducted  with  the  10-deg  and 
1 5-deg  boattail  configurations,  it  was  not  true  for  most  of  the  experiments  with  the  25-deg 
boattail  configuration,  which  had  an  extremely  sharp  shoulder.  j 

If  the  flow  is  "well  behaved,"  as  in  the  former  cases,  the  two  separate  flows  and 
their  mixing  region  are  usually  described  as  indicated  in  Fig.  32.  The  flow  through  the 
nozzle  is  divided  into  an  inviscid  core  region  (A)  and  a thin  boundary  layer  (C),  whereas 
the  external  flow  is  described  in  an  analogous  manner,  the  corresponding  parts  being  (B) 
and  (D). 

Because  of  the  length  of  the  cylindrical  portion  of  the  model  preceding  the  nozzle 
boattail  pictured  in  Fig.  32,  and  because  the  external  flow  is  usually  tripped  somewhere 
near  the  model  nose,  the  boundary  layer  of  the  external  flow  (D)  is  considered  to  be 
turbulent.  This  assumption  is  also  made  for  the  inner  boundary  layer  (C).  Thus,  region 
(E)  is  an  area  of  turbulent  mixing. 

This  relatively  simple  and  straightforward  model  does  not  hold  under  some  conditions, 
particularly  those  situations  that  produce  separated  flow.  Large  closure  rates  of  the  nozzle 
boattail  will  generate  adverse  pressure  gradients  which  can  easily  cause  flow  separation. 

An  increase  in  Mach  number  will  aggravate  the  situation,  particularly  if  it  is  sufficient 
to  produce  local  shock  waves.  Such  flows  are  common  (particularly  with  the  25-deg 
boattail)  but  cannot  be  solved  analytically  at  this  time. 
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5.1.2  Analytical  Model 

The  component  flow  as  described  in  the  previous  section  was  modeled  as  closely 
as  possible.  Although  each  family  of  calculations  will  be  covered  in  detail,  it  may  be 
helpful'  here  to  describe  the  general  solution  procedure. 

Because  a great  deal  of  the  flow  is  subsonic  in  many  cases,  it  would  be  preferable 
to  treat  the  entire  region  of  interest  with  the  Navier-Stokes  equations,  thereby  guaranteeing 
that  the  upstream  influence  of  the  plume  would  be  included  throughout  the  subsonic 
areas.  This  form  of  general  solution  is,  however,  not  yet  available  for  compressible  turbulent 
flows,  and  thus  it  was  decided  to  use  the  component  flow  model  as  shown  in  Fig.  33. 
It  will  be  noted  that  the  only  difference  in  this  model  from  that  in  Fig.  32  is  the  division 
of  the  mixing  layer  into  two  parts  separated  by  the  inviscid  plume  boundary. 

5.1 .3  Computational  Procedure 

The  actual  interaction  between  the  plume  and  the  external  flow  is  so  strong  that 
it  was  deemed  necessary  to  make  the  solution  an  iterative  one,  following  the  outline  in 
Fig.  34.  The  procedure  consists  of  an  inviscid/inviscid  (I/I)  loop  and  a viscous/inviscid 
(V/I)  loop,  in  essence  each  loop  being  an  iteration  within  itself.  The  I/I  loop  is  primarily 
concerned  with  the  consistent  solution  between  the  internal  inviscid  plume  flow  and  the 
external  inviscid  flow,  (A)  and  (B)  in  Fig.  32.  Any  difference  in  pressure  between  the 
external  inviscid  solution  and  the  previously  specified  plume  boundary  pressure  requires 
recalculation  first  of  the  plume  and  then  of  the  external  inviscid  flow  using  the  new  plume 
shape.  In  the  I/I  loop  the  displacement  thicknesses  in  the  boundary  and  shear  layers, 
which  usually  exist  in  all  but  the  first  entry  into  this  loop,  are  assumed  constant  and 
merely  translate  as  the  plume  shape  changes.  When'  reasonable  agreement  is  obtained 
between  the  two  solutions,  the  V/I  loop  is  entered. 

The  V/I  loop  is  concerned  with  adjusting  the  displacement  thicknesses  on  the  nozzle 
boattail  (D)  and  inviscid  plume  surface  (El)  according  to  their  compatibility  with  the 
external  inviscid  flow  (B).  The  inviscid  plume  shape  is  assumed  constant  within  this  loop. 
The  viscous  analysis  which  is  covered  in  Section  5.2.3  and  Appendix  A,  provides  the 
external  displacement  thicknesses.  Any  effect  of  the  internal  boundary  layer  on  the 
afterbody  flow  is  assumed  to  be  sufficiently  small  to  be  ignored.  Upon  leaving  this  loop, 
further  convergence  may  be  obtained  by  reentering  the  I/I  loop.  In  practice  it  has  been 
found  that  this  is  rarely  necessary. 

Inherent  in  the  procedures  described  above  are  the  assumptions  thai  the  boundary 
layers  are  thin  and  that  the  supporting  surfaces  have  little  longitudinal  curvature.  These 
assumptions  may  not  be  strictly  true  for  some  of  the  data  reported  herein,  but  they  are 
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for  the  most  part  valid  for  the  configurations  analyzed,  as  the  flows  remain  attached. 


Further  extension  of  this  procedure  is  expected  into  the  realm  of  separated  flows,  as 
will  be  further  explained  in  Section  5.2.4..  Viscous  Flows. 

5.2  COMPUTATIONAL  TECHNIQUES 

A number  of  different  techniques  were  available  for  the  inviscid  and  viscous  flows, 
and  during  the  initial  computations,  most  of  the  programs  were  used  and  the  results 
compared.  The  strengths  and  weaknesses  of  each  program  gradually  became  apparent,  and 
some  programs  were  eliminated  in  the  subsequent  iterations.  This  was  true  in  particular 
for  the  programs  used  in  the  analysis  of  the  initial  boundary  layers  and  the  inviscid  flow. 
Pertinent  details  of  the  four  main  computations  follow 

5.2.1  Initial  Boundary  Layers 

The  viscous  calculations  in  the  V/I  iterations  were  performed  from  station  130  of 
the  model,  a point  approximately  3.4  radii  ahead  of  the  exit  plane  and  one  radius  ahead 
of  the  initial  point  of  curvature  of  the  nozzle  boattail.  Although  it  was  realized  that  the 
flow  upstream  as  far  as  the  nose  of  the  model  would  be  affected  by  what  was  happening 
at  the  nozzle  exit,  this  station  was  assumed  to  be  sufficiently  upstream  to  be  relatively 
unaffected  when  compared  to  the  variation  near  the  nozzle  exit.  Subsequently,  during 
the  lurse  of  the  iterations,  it  was  found  that  the  displacement  thicknesses  and  surface 
pressures  calculated  in  the  viscous  and  inviscid  portions,  respectively,  of  rhe  iterations 
changed  very  little  during  the  course  of  the  iterations  in  spite  of  large  variations  of  these 
quantities  near  the  nozzle  exit  plane.  Thus,  it  was  felt  that  the  choice  of  MS  130  was 
justified,  although  this  would  not  necessarily  be  the  case  for  the  other  boattail 
configurations. 

In  order  to  determine  the  boundary  layer  at  that  point,  two  boundary  layer  programs. 
Refs.  4 and  5,  were  used  to  calculate  the  boundary  layer  profiles  at  MS  115.46.  The 
two  computations  agreed  well  with  each  other  and  the  wind  tunnel  data,  Fig.  35,  and 
they  were  extended  to  MS  130. 

5.2.2  Plume 

The  plumes  were  calculated  using  the  Lockheed  Method  of  Characteristics  (MOC) 
Program,  Refs.  6 through  9.  The  manner  in  which  this  program  was  used  was  to  specify 
the  Mach  number  distribution  across  the  nozzle,  sonic  in  this  case,  and  the  pressure 
distribution  downstream  from  the  nozzle,  as  well  as  stagnation  properties- of  the  internal 
nozzle  flow.  Calculations  performed  using  the  MOC  program  furnished  the  axisymmetric 
plume  boundary,  hereafter  referred  to  as  the  inviscid  reference  line  (JRL).  This  line  will 
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become  important  in  the  discussion  of  the  viscous  flow  calculations,  Section  5.2.4. 
Characteristic  plume  shapes  for  nozzle  pressure  ratios  of  three  and  seven  are  shown  in 
Fig.  36.  For  the  purposes  of  the  analysis,  the  plume  was  assumed  to  extend  cylindrically 
past  its  point  of  maximum  radius,  as  shown  in  the  figure. 

5.2.3  Inviscid  Flow 

Three  digital  computer  programs  were  used  initially  in  the  inviscid  flow  calculations. 
The  first  program  was  the  Pratt  and  Whitney  Transonic  Program,  Ref.  10,  which  was 
available  and  had  been  used  previously  for  similar  analyses.  Reference  1 1 describes  the 
second  program  used,  a finite-difference  solution  of  the  full  transonic  potential  equation 
by  South  and  Jameson.  This  program  was  modified  slightly  for  our  purposes  so  that  the 
body  shape  could  be  input  directly,  a change  which  gave  the  opportunity  to  then  compute 
the  body  derivatives  using  a number  of  polynomial  or  spline  fits.  The  Douglas-Neumann 
Potential  Flow  Program  was  also  used,  even  though  it  is  primarily  a subsonic  calculation. 
It  does,  however,  incorporate  a Goethert  similarity  transformation  for  high  subsonic  flows. 

In  order  to  compare  the  accuracy,  efficiency,  and  reliability  of  these  three  programs, 
a typical  well-behaved  case  was  chosen  for  the  viscous/inviscid  iterations,  the  15-deg  boattail 
at  = 0.6,  Re/fi  = 2.5  x 106  ft'1  and  NPR  = 3.  The  iterations  were  run  independently 

with  the  three  programs  using  the  same  calculation  procedure  for  the  turbulent 
boundary-lay'r  calculations.  A comparison  of  the  three  "converged"  solutions  is  shown 
in  Fig.  37a. 

The  iteration  technique  used  for  all  the  iterations  except  the  first,  when  the  inviscid 
flow  over  the  body  and  plume  alone  is  calculated,  was  identical  for  all  three  programs 
The  pressure  distribution  from  the  previous  iteration  was  input  into  the  boundary-layer 
program,  and  the  calculated  displacement  thicknesses  were  added  to  the  body.  The  plume 
was  treated  as  an  effective  body  using  a method  which  will  be  described  in  the  next 
section.  This  new  body  shape  was  then  input  to  the  inviscid  program.  Any  discrepancy 
between  the  resulting  pressure  distribution  and  the  previous  one  was  averaged  and  used 
for  the  starting  distribution  for  the  next  iteration.  This  technique  was  used  for  subsequent 
calculations  at  different  Mach  numbers.  All  three  programs  converged  within  four  or  five 
iterations. 

The  general  agreement  among  the  three  programs  was  good,  as  shown  in  Fig.  37a. 
The  disagreement  at  the  afterbody/plume  junction  and  the  Pratt  and  Whitney  Program's 
tendency  not  to  return  to  recovery  pressure  downstream  are  apparent.  Because  of  this 
and  because  true  transonic  capability  was  desired  in  the  inviscid  analysis,  all  subsequent 
iterations  were  performed  with  the  South-Jameson  program,  with  occasional  checks  using 
the  Douglas-Neumann  program. 
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5.2.4  Viscous  Flow 

The  displacement  thickness  over  the  boattail  and,  at  first,  over  the  plume  (considered 
as  a solid  body)  was  calculated  with  the  Kuhn-Nielsen  Turbulent  Separated  Boundary  Layer 
Program  (Ref.  3).  The  program  is  an  integral  technique  utilizing  a 
law-of-the-wall/law-ofthe-wake  velocity  profile  and  an  eddy  viscosity  turbulent  relationship. 
Its  uniqueness  lies  in  the  fact  that  the  mathematical  singularity  occurring  in  the  approach 
to  a point  of  flow  separation  is  bypassed  by  rearrangement  of  the  governing  ordinary 
differential  equations.  Acquired  for  its  separated  flow  capabilities,  the  program  was  not 
used  in  this  mode  for  the  present  analysis,  however,  as  the  1 5-deg  boattail  configuration 
and  plume  analyses  resulted  in  flows  which  did  not  indicate  separation,  as  determined 
by  inspection  of  schlieren  photographs. 

The  problem  of  how  to  handle  the  plume  realistically,  with  both  its  displacement 
and  entrainment  characteristics,  was  accommodated  somewhat  with  the  mixing  analysis 
presented  in  Appendix  A.  This  analysis  showed  that  a technically  valid  technique  was 
to  take  the  value  of  the  displacement  thickness  calculated  at  the  boattail/plume  junction 
and  add  it  to  the  inviscid  reference  line  along  its  entire  length.  Although  this  method 
was  not  entirely  satisfactory,  it  did  eliminate  the  need  to  fair  in  a line  of  arbitrary  shape 
from  the  afterbody  boundary  layer  to  some  point  on  the  inviscid  reference  line,  or  to 
assume  the  plume  to  be  a solid  body.  Further  work  is  continuing  in  this  area  to  include 
the  effects  of  longitudinal  pressure  gradients  and  transverse  (y-direction)  velocities. 

5.3  RESULTS  OF  THE  ANALYSIS 

Figure  37  represents  the  analyses  performed  on  the  1 5-deg  boattail  configuration, 
Re/£  = 2.5  x 106  ft*1,  NPR  = 3,  at  Mach  numbers  of  0.6,  0.8,  and  0.9.  As  pointed 
out  previously  in  Section  5.2.3,  agreement  between  the  three  methods  used  for  the  inviscid 
calculations  was  good  for  the  = 0.6  case,  Fig.  37a.  It  is  apparent  that  overcompression 
is  predicted  near  the  end  of  the  boattail,  a characteristic  of  the  flow  calculation  which 
would  tend  to  make  drag  predictions  highly  optimistic. 

Figures  37b  and  c,  at  Mach  numbers  0.8  and  0.9,  respectively,  show  a continuation 
of  the  ovcrcompression  characteristic.  Although  neither  inviscid  program  accurately  predicts 
the  highly  negative  pressure  at  the  shoulder  of  the  boattail,  the  subsonic  (rather  than 
transonic)  nature  of  the  Douglas-Neumann  program  is  apparent  in  the  solution  at  Mm  = 
0.9  (Fig.  37c). 
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6.0  CONCLUDING  REMARKS 

An  experimental  program  was  conducted  to  investigate  the  interaction  effects  which 
occur  between  the  nozzle  exhaust  flow  and  the  external  flow  field  associated  with  isolated 
nozzle/afterbody  configurations  at  Mach  numbers  0.6,  0.8,  0.9,  0.95,  1.1,  and  1.5. 
Configuration  variables  included  nozzle  boattail  geometries  with  boattail  angles  of  10,  15, 
and  25  deg.  Cold  and  hot  jet  test  techniques  were  used  to  both  simulate  and  duplicate 
the  nozzle  exhaust  flow  for  a sonic  jet  installation.  High  pressure  air  was  used  to  simulate 
the  nozzle  exhaust,  and  an  air/ethylene  combustor  provided  the  means  to  duplicate  the 
nozzle  exhaust.  Nozzle  exhaust  temperature  effects  were  evaluated  over  the  temperature 
range  from  540°R  to  approximately  2,900°R. 

The  significant  results  obtained  from  the  experimental  program  pertain  to  those  effects 
on  boattail  pressure  drag  coefficient  caused  by  flow  asymmetry  (model  support  strut 
induced),  nozzle  base  area,  and  exhaust  plume  temperature.  These  results  are  as  follows: 

1.  Flow  asymmetry  effects  were  Mach  number  and  nozzle  pressure  ratio 
dependent  and  increased  in  severity  as  the  boattail  angle  was  increased. 

For  example,  flow  asymmetry  effects  on  the  25-deg  boattail  configuration 
at  Mach  number  0.95  were  equivalent  to  approximately  16  percent  of  the 
boattail  drag  level. 

2.  Nozzle  base  area  effects  were  significant  and  dependent  on  Mach  number, 
nozzle  pressure  ratios,  and  boattail  configuration.  Increasing  the  nozzle  base 
area  from  0.93  to  6.5  percent  of  the  model  maximum  cross-sectional  area 
resulted  in  increasing  the  boattail  pressure  drag  as  much  as  57  percent  for 
the  15-deg  boattail  at  Mach  number  0.6  and  nozzle  pressure  ratio  3 and 
decreasing  the  boattail  pressure  drag  as  much  as  1 2 percent  for  the  25-deg 
boattail  at  Mach  number  0.95  and  nozzle  pressure  ratio  3. 

3.  The  differences  obtained  in  boattail  pressure  drag  between  the  cold  jet 
simulation  and  hot  jet  duplication  results  were  significant  at  nozzle  pressure 
ratios  representative  of  turbofan  and  turbojet  engines  at  subsonic  Mach 
numbers.  Adjusting  the  cold  jet  nozzle  pressure  ratio  to  correct  for  the 
changes  in  the  exhaust  plume  specific  heat  ratio  with  temperature  did  not 
account  for  the  differences  observed.  At  Mach  number  0.8  and  for  a nozzle 
pressure  ratio  of  3,  adjusting  the  15-deg  boattail,  cold  jet  results  to  a result 
representative  of  an  engine  operating  with  a plume  temperature  of  2,900°R 
would  result  in  a 45-percent  error  (increase  in  drag)  in  boattail  drag. 
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In  a continuing  effort  to  improve  nozzle  afterbody-related  theoretical  prediction 
techniques,  the  experimentally  measured  pressure  distribution  (top  row  of  pressure  orifices) 
on  the  15-deg  boattail  configuration  at  Mach  numbers  0.6,  0.8,  and  0.9  were  compared 
with  theoretical  results  obtained  from  the  South-Jameson,  Douglas-Neumann,  and  Pratt 
and  Whitney  inviscid  programs  in  conjunction  with  the  Kuhn-Nielson  separated 
boundary-layer  program.  The  analyses  were  limited  to  that  for  a cold  jet  operating  at 
a nozzle  total  pressure  ratio  of  3 and  were  performed  using  a viscous/inviscid  iterative 
procedure.  All  three  computer  programs  predicted  more  recompression  near  the 
boattail/plume  junction  than  was  actually  measured.  Further  work  is  continuing  to  remedy 
this  problem,  particularly  through  more  sophisticated  mixing  analyses.  The  programs  predict 
with  reasonable  accuracy  the  low  pressures  encountered  at  the  boattail  shoulder;  the 
South-Jameson  program  was  the  preferred  program  since  it  predicted  results  nearer  the 
actual  measurements.  Predictions  of  pressure  level  at  other  points  along  the  boattail  agreed 
reasonably  well  with  the  actual  measurements  with  results  improving  as  the  free-stream 
Mach  number  was  decreased  from  sonic  conditions. 
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Figure  1.  Basic  model  dimensions  and  location  in  test  section. 
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Figure  3.  Concluded.  i 
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Figure  11.  Hot  flow  model  exhaust  total  temperature. 
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Fiaure  13.  Continued. 
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Figure  14.  Continued 
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Figure  15.  Continued. 
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Figure  15.  Continued. 


f..  M„=  1.50 
Figure  15.  Concluded. 
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AH  DC-TR  -76-1 02 


10-deg  boattail,  NPR  = 6.83 


15-deg  boattail,  NPR  = 5.7 


25-deg  boattail  NPR  = 6.99 
a.  M = 0.60 

oo 

Figure  16.  Schlieren  photographs  of  10-,  15-,  and  25-deg  boattail 
geometries  at  selected  Mach  numbers. 
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Variation  of  boattail  pressure  drag  coefficient  with  Mach  number. 
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Figure  19.  Supersonic  Mach  number  effects  on  boattail  pressure  distribution. 
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Figure  20.  Flow  asymmetry  effects  on  10 -deg  boattail  pressure  drag  coefficient. 


Bottom 


asymmetry  affects  on  25*deg  boattail  pressure  dreg  coefficient. 


Continued 


22.  Concluded. 


a.  10-dag  boattail 

mmetry  affects  as  a function  of  Mach  number. 


Continued 


c.  M„  = 1.50 
Figure  24.  Concluded 


asymmetry  effects  on  15-deg  boattail  pressure  distribution 


Concluded 


Concluded 


a.  M.  = 0.60 

Figure  28.  Exhaust  plume  temperature  effects  on  25-deg  boattail  configuration. 
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CDPT 


Sym  TTJ , °R  PN/TP 
A 540  2 001 

B 2,150  4601 

C 3,025  1901 


0 10  20  30  40 

DELH,  deg 


d.  M„  = 0.95 
Figure  29.  Continued 


AEDC-TR-76-102 


0.22 

CDPT 

0.20 

0.18 

0.16 

0.14 

0.12 


d.  - 0.95 
Figure  30.  Continued 
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S^m 

Boat  tail,  deg 

ABASE/A 

A 

15 

0.0093 

B 

25 

0.0093 

0 

15 

0.0650 

□ 

25 

0.0650 

S^m 

Boat  bail,  dee 

ABASE/A 

A 

15 

0.0093 

B 

25 

0.0093 

0 

15 

0.0650 

□ 

25 

0.0650 

Schematic  of  the  various  flow  regions  involved  in  the  afterbody  problem 


Figure  33.  Schematic  of  the  various  flow  region*  involved  in  the  modeling  of  the  afterbody  problem. 
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Figure  34,  General  iteration  flow  plan. 
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Table  1.  External  Pressure  Orifice  Locations 
a.  10-deg  Boattail  Configuration 


X 

X/D 

R 

Pr 

sssure  Orifices 

Top 

Bottom 

Side 

-15. 581 

ESS 

4.930 

■M 

607 

— 

-11.211 

WSEm 

608 

mm 

- 6.831 

604 

609 

BBfl 

- 2.461 

RSH 

605 

610 

--- 

1.  522 

0.154 

1 

* 

401 

442 

464 

3.914 

0.397 

4.799 

402 

443 

465 

4.874 

0.494 

4.706 

403 

444 

466 

5.683 

0.576 

4.613 

404 

445 

467 

6.  509 

0.660 

4.514 

405 

446 

468 

7.178 

0.728 

4.426 

406 

447 

469 

7.893 

0.801 

4.332 

407 

448 

470 

8.570 

0.869 

4.229 

408 

449 

471 

9.292 

0.942 

4.120 

409 

450 

472 

10.023 

1.017 

4.000 

410 

451 

473 

10.506 

1.075 

3.896 

411 

452 

474 

11.186 

1.134 

3.780 

412 

453 

475 

11.789 

1.196 

3.646 

413 

454 

476 

12.336 

1.251 

3.513 

414 

455 

477 

12.834 

1.302 

3.379 

415 

456 

478 

13.283 

1.348 

3.246 

416 

457 

479 

13.750 

1.395 

3.094 

417 

458 

480 

14.180 

1.438 

2.941 

418 

459 

481 

14.628 

1.484 

2.773 

419 

460 

482 

15. 105 

1.532 

2.584 

420 

461 

483 

15.670 

1.589 

2.342 

421 

462 

484 

16.210 



1.644 

2.134 

422 

463 

485 

X - Location  Forward  and  Aft 
MS  130.471 

D - Model  Maximum  Diameter,  9.  86  in. 
H - Radius 


View 

Looking 

Upstream 
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Table  1.  Continued 
b.  15-deg  Boattail  Configuration 


Pressure  Orifices 


•15. 581 
•11.211 

- 6.831 

- 2.461 

1.  730 

2.  725 

3.  727 

4.725 

5.725 

6.725 

7.  503 

8.  156 

8.  707 
9.302 

9.  8 CO 
10.294 
10.764 

11. 154 
11. 565 
11.958 
12.385 
12. 787 

13. 155 
13.657 
13.960 
14.581 
15.069 
15. 560 
16. 120 


•1. 580 
•1.137 
•0.693 
-0.250 
0.  175 
0.276 
0.378 
0.479 
0.  581 
0.682 
0.761 
0.  827 
0.883 
0.943 

0. 999 
1.044 
1.092 

1.  131 
1.173 
1.213 
1.256 
1.297 
1.334 
1.385 
1.416 
1.479 
1.  528 
1.  578 
1.635 


rt 

Top 

Bottom 

Side 

4.9 

po 

602 

607 

- 

603 

608 

— 

604 

609 

— 

605 

610 

— 

1 

f 

401 

— 

— 

4.925 

402 

— 

— 

4.921 

403 

— 

— 

4.920 

404 

— 

— 

4.919 

405 

— 

— 

4.884 

406 

442 

464 

4.795 

407 

443 

465 

4.693 

408 

444 

466 

4.  591 

409 

445 

467 

4.470 

410 

446 

468 

4.352 

-11 

447 

469 

4.245 

412 

448 

470 

4.  117 

413 

449 

471 

4.004 

414 

450 

472 

3.  882 

415 

451 

473 

3.756 

416 

>52 

474 

3.615 

417 

3 

475 

3.447 

418 

45-. 

476 

3.345 

419 

455 

477 

3. 160 

420 

456 

478 

3.042 

421 

457 

479 

2.795 

422 

458 

480 

2.  584 

423 

459 

481 

2.364 

424 

460 

482 

2.  144 

425 

461 

483 

X - Location  Forward  and  Aft 
MS  130.471 

D - Model  Maximum  Diameter, 
9. 86  in. 

R - Radius 


Side  -| 


View 

Looking 

Upstream 


Bottom 
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Table  1.  Concluded 
c.  25-deg  Boattail  Configuration 


X 

X/D 

R 

Pressure  Orifices 

Top 

Bottom 

Side 

-15. 581 

-1. 580 

4.930 

602 

607 

-11.211 

-1. 137 

603 

608 

--- 

- 6.831 

-0. 693 

604 

609 

— 

- 2.461 

-0.250 

605 

610 

— 

2.433 

0.247 

401 

— 

— 

3.433 

0.348 

402 

— 

— 

4.433 

0.450 

403 

— 

— 

5.433 

0.  551 

404 

— 

— 

6.433 

0.  652 

405 

— 

7.433 

0.754 

406 

— 

— 

8.433 

0.855 

407 

— 

— 

9.431 

0.956 

’ 

f 

408 

— - 

— 

10.430 

1.058 

4.886 

409 

442 

464 

10.730 

1.088 

4.800 

410 

443 

465 

11.040 

1. 120 

4.  672 

411 

444 

466 

11.358 

1.152 

4.509 

412 

445 

467 

11.682 

1.185 

4.342 

413 

446 

468 

12.029 

1.220 

4.  167 

414 

447 

469 

12.381 

1.256 

3.986 

415 

448 

470 

12.761 

1.294 

3.798 

416 

449 

471 

13.151 

1.334 

3.602 

417 

450 

472 

13. 562 

1.375 

3.403 

418 

451 

473 

13.995 

1.419 

3.200 

419 

452 

474 

14.461 

1.467 

2.982 

420 

453 

475 

14.969 

1.  518 

2.  747 

421 

454 

476 

15.516 

1.574 

2.485 

422 

455 

477 

16.  138 

1.637 

2- 

198 

423 

456 

478 

X - Location  Forward  and  Aft 
MS  130.471 

D - Model  Maximum  Diameter, 
9.86  in. 

E - Radius 


Top 


Bottom 


View 

Looking 

Upstream 
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Table  2.  Cold  Flow  Nozzle  Assembly  Instrumentation 


Total  Pressure  Hake  (MS  133.  587) 


Pressure  Position  (Radius,  in.) 


PTJ  1 
2 

3 

4 

5 

6 

7 

8 
9 

" 10 


Nozzle  Centerline 
0,  389 

0. 778 

1.  167 
1.  556 

1.  945 
2,334 

2.  723 

3.  112 
3.  501 


Wall  Static  Pressures  (MS  133.  57) 


Pressure  Position,  deg 


PW  1 
2 
3 

" 4 


45 

135 

225 

315 


0 


View 

(+)  Looking 
Upstream 


Total  Temperature  Probes  (MS  135,70) 


Temperature 

TD  1 
TD  2 


Position,  deg 

120 

240 
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Table  3.  Hot  Flow  Nozzle  Assembly  Instrumentation 


Chamber  Wall  Static  Pressures  (MS  136. 47) 


Pressure 


PW  1 
2 

3 

4 


Position,  deg 

45 

135 

225 

315 


o 


View 

Looking 

Upstream 


Nozzle  Exit  Static  Pressures  (MS  146.  895) 


Pressure  Position,  deg 

PE  1 45 

2 135 

3 225 

11  4 315 


Chamber  Liner  Temperatures  (MS  146.  895) 

Temperature  Position,  deg 

TL  1 0 

TL  2 180 


I 


1 

i 

135 

> ' 
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APPENDIX  A 

PLUME  DISPLACEMENT  THICKNESS  CALCULATION 


Assume  that  the  static  pressure  is  everywhere  constant  and  that  the  flow  direction 
is  parallel  to  the  inviscid  reference  line  (IRL);  that  is,  transverse  (y-direction)  flow  may  be 
neglected. 


Mass  is  conserved  above  the  DSL 

y„-y 


U ] SO  * U 

f p u dy , = / p u dy 

° yp 

5 1 y «•->•»  yu 

/ Pudyj  + / PjUjdy,  = f P udy 

° SI  XD 


(1) 


/ + P]«)(yu  - >'«  - sp  = J pudy 

0 yp 

Apply  the  definition  of  Sv  at  the  downstream  position  assuming  the  DSL  as  the  surface, 
as  follows: 

* y u 


- / "(i  - 

>p  \ Pl"V 


/ y = - p i«l(«p  - y„  + y1}) 


(2) 
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Combining  (1)  and  (2)  gives 

/ pudy,  + plUj(yD  + S - y^,  - 5,)  = 0 
0 

Apply  the  definition  of  6*  to  the  upstream  position,  as  follows: 


/ pufyi  = PjUi^!  - 5 j) 


O 

Substitution  of  (4)  into  (3)  gives 

S*P  = 8\  + (y-  ~ >’n^ 

The  displacement  relative  to  the  IRL  boundary  is,  by  definition 

5irl  * (>'n  - yJ  + 

Combining  Eqs.  (5)  and  (6)  yields 


5 


♦ 

1 


(3) 

(4) 

(5) 

(6) 


Thus,  the  "plume  displacement  thickness"  relative  to  the  IRL  is  simply  the  displacement 
thickness  of  the  boundary  layer  at  the  afterbody/plume  junction. 
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NOMENCLATURE 

A Model  cross-sectional  area,  0.5303  ft2 

ABASE  Nozzle  annular  base  area,  0.03450  ft2  (hot  flow  nozzle)  and  0.00492  ft2 

(cold  flow  nozzle) 

AT  Nozzle  throat  area,  0.05683  ft2  (hot  flow  nozzle)  and  0.08648  ft2  (cold 

flow  nozzle) 

CDND  Nozzle  discharge  coefficient 

CDPX  Boattail  pressure  drag  coefficient,  X = T(top),  B(bottom),  or  S(side) 

CPj  Pressure  coefficient,  (P;  - P^)/qM 

CP(l.O-X)  Pressure  coefficient  for  10-deg  boattail  configuration,  X = T(top), 

B(bottom),  or  S(side) 

CP(15-X)  Pressure  coefficient  for  15-deg  boattail  configuration,  X = T(top), 

B(bottom),  or  S(side) 

CP(25-X)  Pressure  coefficient  for  25-deg  boattail  configuration,  X = T(top), 

B(bottom),  or  S(side) 

CP*  Pressure  coefficient  at  sonic  Mach  number 

D Model  maximum  diameter  0.82166  ft 

DELH  Nozzle  jet  initial  inclination  angle,  deg 

DSL  Dividing  streamline 

FTA  Fuel/air  ratio 

IRL  Inviscid  reference  line 

I/I  Inviscid/inviscid  iteration  loop 

K-y  Specific  heat  ratio  function,  0.77658  Oy)0-5 exp 

Ma  Air  mass  flow,  Ibm/sec 

Mf  Fuel  mass  flow,  Ibm/sec 
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MOC  Method  of  characteristics 

MS  Model  station,  in. 

Mm  Free-stream  Mach  number 

NPR  Nozzle  total  pressure  to  free-stream  static  pressure  ratio 

NPRB  Nozzle  total  pressure  to  free-stream  static  pressure  ratio,  backup 

Pi  Static  pressure,  psfa 

Pw all  Nozzle  flow  duct  static  pressure  forward  of  throat  contraction,  psfa 

P^  Free-stream  static  pressure,  psfa 

PN/TP  Part  number/test  point 

PTave  Nozzle  flow  duct  rake  average  total  pressure,  psfa 

PTJ  Nozzle  flow  duct  calculated  total  pressure,  psfa 

PTX  Nozzle  flow  duct  rake  total  pressure,  psfa 

(PT/P)i  Isentropic  pressure  ratio 

Free-stream  dynamic  pressure,  psf 
Re/fi  Free-stream  Reynolds  number  per  foot,  1/ft 

RET  Nozzle  throat  Reynolds  number 

R Cold  flow  nozzle  duct  radius  0.3004  ft 

R,  Nozzle  boattail  radius,  ft 

TTJ  Nozzle  exhaust  average  total  temperature,  °R 

Ue  External  velocity 

U,  Local  velocity,  ft/sec 

U„  Free-stream  velocity,  ft/sec 

V/I  Viscous/inviscid  iteration  loop 
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X/D  External  pressure  orifice  location  relative  to  MS  130.471 

Yj  Transverse  coordinate 

Y/R  Total  pressure  probe  location  from  nozzle  flow  duct  wall 

7 Ratio  of  specific  heats 

5 Boundary-layer  thickness,  ft 

5*  Boundary-layer  displacement  thickness,  ft 

p Density,  lbm/ft3 


